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Abstract 



We study here the behaviour of non-spacelike geodesies in dust collapse models 
in order to understand the casual structure of the spacetime. The geodesic families 
coming out, when the singularity is naked, corresponding to different initial data are 
worked out and analyzed. We also bring out the similarity of the limiting behaviour 
for different types of geodesies in the limit of approach to the singularity. 

1 Introduction 

The gravitational collapse of a spherically symmetric dust cloud is described by the Tolman- 
Bondi-Lemaitre (TBL)[|TJ models. These models have been extensively studied for the va- 
lidity of the cosmic censorship conjecture @]. In particular, it is known || now that 
depending upon the initial conditions, which are denned in terms of the initial density and 
velocity profiles from which the collapse develops, the central shell-focusing singularity at 
r = can be either a black hole, or a locally or globally naked singularity. Further, it has 
been shown recently that this singularity is always gravitationally strong || along a family of 
timelike geodesies, indicating the physical nature and significance of the singularity. There 
have also been recent studies which examine the stability and other aspects in this case 0. 

One may thus state that the final fate of gravitational collapse is reasonably well- 
understood for the dust collapse models. The same however cannot be said here about 
the causal structure of the singularity, on which several aspects are still not clear, especially 
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when it is naked. It is the behaviour of null geodesies families, closer to the singularity, that 
give us an idea of the causal structure there. Most of the earlier works on dust collapse have 
concentrated mainly on establishing the occurrence of black holes and naked singularities in 
TBL models under various sets of initial conditions. For the case of a black hole forming in 
dust collapse, the behaviour of null geodesies families around and close to the event horizon 
is rather well-known and fully explored. However, this is not the case for the behaviour of 
null geodesies in the vicinity of a naked singularity, which is largely unknown. We take up 
here a study of the behaviour of radial geodesies in the neighbourhood of a naked singularity. 
We bring out the interesting behaviour of null trajectories near naked singularities, and it is 
in this sense that we explore this causal structure here. Understanding such issues is clearly 
going to be important if there is any interesting physics to come out due to the existence of 
a naked singularity. Even from the perspective of obtaining a suitable formulation of cosmic 
censorship, such an understanding should be helpful and necessary. Our purpose here is to 
investigate the nature and causal structure of this singularity by means of examining the 
behaviour of null as well as timelike geodesies families in the TBL models. This provides a 
better understanding of the structure of this singularity, and clarifies several related issues. 

It is seen from our considerations here that there is a single ingoing radial null geodesic 
(RNG) terminating at the singularity with a well-defined tangent, and that there is a single 
RNG (the Cauchy horizon) coming out along one direction, and a family of RNGs coming 
out in the other direction in all the cases in which the singularity is naked. We also show 
that there is a family of radial timelike and radial spacelike geodesies coming out of the 



3 



singularity with well-defined tangents. We explicitly show these results for cases which were 
not studied so far. In some other cases, where one had some idea on the behaviour of the 
RNGs, the earlier results can again be reconfirmed from the present calculations. 

In particular, with our scaling when the physical radius (R) along the geodesies is pro- 
portional to a power of the comoving radius (r) less than 3, it is seen that there is only one 
ingoing radial null geodesic (RNG) terminating at the singularity, and one RNG coming out 
along the direction of Cauchy horizon as stated above. Further, there is a family of radial 
null and timelike geodesies coming out of the singularity along the direction of the apparent 
horizon. In fact, earlier it was claimed in |4]], giving an evidence for the numerical results 
obtained there, that in certain cases there is a family of RNGs coming out of the singularity 
along the apparent horizon. Here we explicitly bring out the existence of such a family, and 
show that the results can be generalized for timelike and spacelike geodesies as well. 

The plan of the paper is as follows. Section 2 introduces the TBL model and the gravita- 
tional collapse. In Section 3 we study the various singular geodesies, and geodesic families. 
In the concluding section 4 we discuss the overall scenario, and the possibility of generalizing 
these results. 
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2 The TBL model and gravitational collapse 



The metric for the TBL models in comoving coordinates is written as, 



ds * = _ dt 2 + ^jdr 2 + R 2 dn 2 . (1) 



The energy-momentum tensor is that of dust, 



T ij = e5t5i, c = c(«,r) = -^, (2) 



where e is the energy density, and the area radius R = R(t, r) is given by 



R 2 = f(r) + (3) 



Here the dot and prime denote partial derivatives with respect to the coordinates t and r 
respectively, and for the case of collapse we have R < 0. The functions F and / are called 
the mass and energy functions respectively, and they are related to the initial mass and 
velocity distribution in the cloud. 

For a marginally bound cloud (/ = 0), the integration of equation (3) gives 



' IJr} w (4) 
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where £o( r ) is a constant of integration. Using the coordinate freedom for scaling the radial 
coordinate r we set, 



R(0,r) = r, 



(5) 



which gives, 



to(r) 



2r 3/2 



(6) 



3^ 



At the time t = t (r) the shell labelled by the coordinate radius r becomes singular because 
the area radius R of the shell becomes zero then. Here we consider only the situation where 
there are no shell-crossings in the spacetime. A sufficient condition for this is that the 
density be a decreasing function of r, which may be considered to be a physically reasonable 
requirement, because for any realistic density profile the density should be higher at the 
center, decreasing away from the center. The ranges of coordinates are given by, 



The quantity R', which is also needed later in the equation of RNGs to check the visibility 
or otherwise of the central singularity, can be written as, 



Though we consider here the marginally bound case / = for the sake of clarity and 
simplicity, the results are easily generalized for non-marginal case. 



< r < oo, — oo < t < to(r). 



(7) 




(8) 
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3 Non-spacelike geodesies in TBL models 



In this Section we will first study the basic equations for geodesies and write down the root 
equation for existence of geodesies coming out of the singularity. Then we show the existence 
of a family of RNGs along the apparent horizon direction, and finally we study the timelike 
geodesies in some detail. 

3.1 Basic geodesic equations 

The equations for the radial geodesies in TBL model can be written as, 



where K l is the tangent vector to the geodesic, k is affine parameter. In our notation if there 
are two signs upper sign always represents the equation for outgoing geodesic while the lower 
one represents ingoing geodesies. The function P has to satisfy the differential equation, 



dt 
dk 



P 



(9) 




(10) 



dP 
~dk 



+ (P 2 + BR 2 ) 



YRR' R 2 \ 



=F VP 2 + BR 2 ^- - BtJf/r = 0. 



(11) 



This gives, 



In the above equations B = for null geodesies, B = — 1 for timelike geodesies, and B = 1 for 
spacelike geodesies. We have dt/dr > for the outgoing geodesies, and they are represented 
by the upper sign in above equations. We introduce new variables u = r a (where a > 1 
is a constant), X = R/u, and write all the quantities in terms of these variables. Then, in 
the limit of approach to the singularity, using the l'Hospital's rule we can write the limiting 
value of X as, 

D JD 

X = lim X= lim -= lim — = U(X ,Q) 

or, 



U(X,0)-X = V(X) = 0. (14) 

where U(X, u) = dR/du (along the geodesies). If the above equation has a real positive root 
X = Xq then the singularity is at least locally naked 0. 

Near the singularity we assume the form of the mass function to be, 

F = F r 3 + F n r 3+n + higher ordered terms. 

The first term on the right hand side in the above equation is required by the condition 
of regularity of the density profile on the initial surface with our scaling. We need to keep 
the first two lowest order terms in the mass function to get the required behaviour of the 
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geodesies near the central singularity. Here we will mainly consider n < 3, i.e. the cases 
with a = 1 + 2n/3 < 3 in the earlier papers ||. In these cases the singularity is always 
naked if the density decreases as we go away from the centre. The value of the (larger) real 
positive root for (|l4f) is given by X = (—F n /2F ) 2 ^ 3 [3], and it is the same for both ingoing 
and outgoing geodesies apart from the two special cases discussed in [|J. 

It can be shown (see e.g. Joshi and Dwivedi in Ref. ||; see also f|) that as (dV/dX)\ x =x 
= Hq < 1 there is only one RNG terminating at the singularity with the larger root Xo as a 
tangent. Similarly in these cases it can be shown that there is only one ingoing radial null 
geodesic that terminates at the singularity as the V(X) equation is the same for both the 
ingoing and outgoing radial null geodesies, and the geodesic equations differ only in higher 
ordered terms. We note that here we mainly study the outgoing geodesies. 

3.2 The family of singular geodesies 

Consider now the collapse with initial densities as generated by the mass function as above, 
with a < 3. This corresponds to the initial density profiles as given by either p(r) = po + pir, 
or p(r) = po + P2 r2 , with pi,p2 < 0. We now analyze in these cases the structure of the 
geodesies coming out from the singularity to show that there is such a family of outgoing 
RNGs terminating at the singularity in the past along the apparent horizon. 

Firstly, let us see in a transparent manner how such a behaviour is possible for the null 
geodesies coming out, and terminating in the past at the singularity. For such a purpose, 
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take the equation of the curve to be R(r) = br ai (where ati > 1 and b are constants for which 
the value is to be fixed later) to the lowest order, and we see the conditions for it to be null 
as r — > 0. We assume that F has form F = F Q r 3 + F n r 3+n . Using the TBL solution we get 
along this curve, 



dt nF„ „ , , _ 6 3/ 



2 



dr 3Fn 3/2 



^ r -i_( ai _l)^_ r |K-i)-i 



and 



' V r 



R' = for" 1-1 - n ^ n r^^ 1 -"!). 



3F Vb 

Note that we have kept only the lowest order terms which will be required for the analysis 
near the central singularity for trajectories which terminate at the central singularity. One 
can take care of more general cases by (explicitly) keeping to( r ) i n the above equations. For 
the curves to be singular we need a\ > 1 with our scaling. The conditions for this curve to 
be null is 



dt 
dr 



R' 

When n = 3, we see that this is possible only when a.\ = 3, and all the powers of r in 
the numerator and denominator become equal and we get the earlier roots equation [[jj], i.e. 
V(b) = V(X) = for nakedness. When n < 3 we can have two cases. If a\ = 3, then the 
first term in dt/dr and the second term in R' dominate, they have equal powers of r and 
we get b = F Q in the limit for the outgoing RNGs. This shows that in the limit in the past 
the outgoing null geodesies have a similar behaviour to that of the apparent horizon. The 
second possible case is ai = 1 + 2n/3. In this case the two terms in dt/dr have an equal 
power, and the two terms in R' also have equal powers. But the power of r in dt/dr is less, 
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so we need that the coefficient should vanish and again we get the earlier roots equation 
Xq = (—F n /2Fo) 3 / 2 [?[, and we need higher order correction terms to R(r) to cancel the 
powers of r. It has been shown earlier that there is only one outgoing or ingoing RNG along 
this direction because the value of ho = (dV/dX)\x=x is less than 1. For all other values of 
«i we see that the two terms in R! and the two terms in dt/dr have different powers, and 
also the lowest powers of r in R' and dt/dr are different. So we cannot have singular null 
geodesies with other values of ax possible. Thus we again see that for n < 3 the singularity 
is always at least locally naked. 

Now, we try to check in a simple way when we have a root to our root equation, whether 
a family of geodesies can terminate at the singularity with the given root as a tangent. For 
simplicity and clarity we will discuss only the outgoing geodesies, but with the same method 
one can easily analyze the ingoing geodesies as well. 

Up to the lowest order, for the singular geodesies, the value of the root (tangent) is decided 
by the self-consistency of the differential equations for geodesies; the constant of integration 
(corresponding to the family of geodesies if it comes out along the given root direction) can 
come through only higher ordered (additive) terms. To check for such a family to exist, we 
assume that along the geodesies the area radius R has the following form (here first we will 
consider the normal root direction), 

R& X r a + Dw{r), 
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where D is considered to be the constant of integration which can label different geodesies 
and w(r) is a function through which the behaviour of the family comes out. 



Now with the form of mass function considered above (F = F r 3 + F n r 3+n ), both the 
sides of the geodesic equation (dR/dr) can be written to the lowest order of these terms as 



(with a = 1 + 2n/3 in this case), 

dw(r) 



aXor^ 1 + D 



dr 



/ZiL r W3 + fF° Dw ( r ) r -n 



X 



Xor 2n/3 + 



Dw(r) 



nF n 



r 2n/3 + 



nF n Dw(r) 
QF X 3/2 r 



(15) 



So, basically after satisfying the root equation, the differential equation for w(r) becomes, 



dr 



1 — [T^ r l ~ n /^ ) ( 1 + 



X 



nF n 



6F X, 



:s/2 



+ 2^v5' x ° 



nF n 



3FqVXo 



,1— ji/3 



U dr dX 



(16) 



dr 



dX 



X=X 
X=X y 



We define, 



ho = 1 + 



(17) 



Now let us first consider the n < 3 (i.e. a < 3) case. In this case the root is given by, 
X^ 2 = —F n /2F , and so after canceling the terms involving the root (as X satisfies the 
root equation for the geodesies) we get near the central singularity, 



D 



dw(r) _ n w(r) 
dr 



D- 



nF n 



6F Xo 



3/2 
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i.e. 



dw(r) w(r) 



dr 



(l-n/3). 



This, after integration gives, 



w(r) oc r 1 n/3 . 



But as 1 — n/3 < 1, to(r) goes to zero slower than r, so we have only one RNG (D = 0) 
coming out along this direction. Thus we get the similar result to that obtained by the 
earlier method 131. 



Now let us consider the n = 3, i.e a = 3 case. In this case, when the singularity is naked 
we have two real positive roots to the root equation. After cancelling the terms satisfying 
the root equation and retaining the lowest powers of w(r), we get the equation, 



D 



dw{r) 
dr 



D- 



w r 



nF„ 



3-Fov Xq 



which after integration gives, 



w(r) oc r 



1 1 i 



6F X^ /2 / 2X V X 



(19) 



So, to have a family of geodesies coming out along the tangent direction X we need, 



ho 



Xj\ Q Fo X 3 /2 J 2X VX 



nF„ 



3F Q ^/X Q 



> 1. 



(20) 



This is the same result as shown earlier by Joshi and Dwivedi. In this case when the 
singularity is naked there are two real positive roots for the root equation, V(X) = 0. So 
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along one root h — 1 is negative and along the other it is positive. As the coefficient of 
highest power of X in V(X) is negative, along the larger root h — 1 is negative while along 
the smaller root it is positive. That means along the larger root we have h < 1, and we can 
have only one RNG coming out along that direction (with D — 0). Along the smaller root 
h > 1 and we have an infinite family of RNGs coming out along this direction. 

Now let us see what happens in the n < 3, i.e a < 3 cases, along the apparent horizon 
direction. In this case we need to keep the two lowest order terms to make sure that we have 
the apparent horizon as a tangent near the singularity. In this case, to the lowest required 
order, the equation becomes, 

9 ^dw(r) 9 w(r) 1 / nF„ „ ,\ 

This, after integration gives, 

nFn 1 

W (r)oce«o /2 ^. (21) 

That means w(r) goes to zero exponentially and we have a family of RNGs coming out of 
the singularity with the apparent horizon kind of behaviour. 

The Fig. 1 and Fig. 2 show the plots of the geodesies equations. Fig. 1 shows the plots 
for F = F r 3 + F x r A ", and Fig. 2 shows the plots for F = F r 3 + F 2 r 5 . We have also plotted 
the apparent horizon, and the analytical form of the geodesic equation near the singularity 
with first correction to F [4], i.e. 
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A 0.015 - 



R 0.010 




Apparent horizon 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.3 

r • 

(a) 



R=F 



R=F-18F 7, V- n /(nF n ) 



0.01 0.02 0.03 0.04 0.05 

r • 



(b) 



Figure 1: In this models density goes to zero smoothly at the boundary. The energy and 
mass functions are / = 0, F = r 3 — 2.5r 4 . a) Plots of different outgoing radial null geodesies 
terminating at the singularity in past. We clearly see the R = F kind of behaviour for these 
geodesies near the singularity, b) graph showing error Ractual ~ Ra pp roximate against r near the 

^actual 

singularity for R = F curve and the first corrected curve R = F 



18K? 



nF n 



o r 6— n 



R = F — l*S! r «-» 
nF n 



It matches very well the actual geodesic equations near the singularity. Fig. 1(b) and Fig. 
2(b) show the error in analytical approximation against the numerical solution at different 
values of r near the singularity. We see that the first corrected result to F matches much 
better with the actual trajectories near the singularity. The plots also give strong numerical 
evidence that there is a family of RNGs meeting the singularity along the apparent horizon. 
The deviation of geodesies from each other seems to be very drastic as they move out. In 
other words, near the center all the singular geodesies converge very fast and fall on each 



other as predicted by our analytic calculation. 
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0.14 




Figure 2: Density goes to zero smoothly at the boundary The energy and mass functions 
are / = 0, F = r 3 — 2.5r 5 . a) Plots of different outgoing radial null geodesies terminating 
at the singularity. We clearly see the R = F kind of behaviour for these geodesies near the 
singularity, b) The graph showing error Ractual ~ Ra pp roximate against r near the singularity for 

^actual 

R = F curve and the first corrected curve R = F g— r 6_n . 



3.3 The timelike geodesies 



Coming to timelike geodesies, by using the set of differential equations (0-0) for the geodesies, 
and using self-consistency requirements, the approximate solution for the timelike geodesies 
near the central singularity along the normal root direction (R = X r a ) can be written as, 



K l = — — ± r 1- "/ 3 — — + higher order terms, (22) 

2C 2C |I 



R(t,r) =X r 1+2n / 3 , (23) 
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Here C 7^ 1 is a constant labeling different geodesies corresponding to the particles having 
different energies at a given point in spacetime, and the signs ± represent outgoing and 
ingoing geodesies. The case C = 1 was discussed in [5]. So we see that when C ^ 1 the 
timelike and spacelike geodesies have the same limiting tangent at the singularity as that of 
the RNGs terminating at the singularity. Thus all the geodesies go along the same direction, 
which we would not have expected to happen typically. 

Let us now look for timelike radial geodesies coming out of the central (r = 0) sin- 
gularity along the apparent horizon, if they exist, i.e. with the behaviour R = F r 3 + 
higher ordered terms. Assuming such geodesies to exist, and using the set of differential 
equations (^| - [H]) for the timelike geodesies and solving them in an approximate manner 
near the singularity, we get, 



K l = C Exp 



~ nFn _ r n-3 



3(3 - n)F 5/2 



(24) 



Here different values of C represents different geodesies corresponding to the particles having 
different energies at a given point in the spacetime like in the earlier equation. We see that 
K l and K r blow up exponentially, making K l / K r = 1, and to the lowest orders, near the 
singularity the dR/dr equation becomes the same as that of the null geodesies. Along any 
other directions (i.e. apart from the larger root, which is the Cauchy horizon direction, 
and the apparent horizon direction, and the two special cases discussed earlier 0), we get 
contradiction while solving the set of differential equations for the geodesies. That means 
these are the only possible behaviours for the geodesies near the singularity. And because K l 
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diverges very rapidly near the central singularity, all the arguments and calculations given 
above for null geodesies also hold for timelike geodesies, i.e. families of timelike and null 
geodesies come out of the singularity along the apparent horizon direction near the central 
singularity. If one also solves the set of differential equations for spacelike geodesies in a 
self-consistent manner, then one can see that even the spacelike geodesies have a behaviour 
very similar to that of null geodesies near the central singularity. 

Fig. 3 shows the timelike geodesies. Within a given range, timelike geodesies starting 
from the singularity can meet the boundary of the cloud at the same time with different 
velocities, and there is also a family of such geodesies which meet the boundary of the cloud 
at different times as shown in this Figure. This implies that in some region, at every point 
in the spacetime we have timelike geodesies coming out of the singularity with apparent 
horizon as tangent at the singularity and meeting that point with different velocities. 

Though we have considered here the marginally bound case / = for the sake of clarity 
and simplicity, the results are easily generalized for non-marginal cases as well. They essen- 
tially depend upon the value of the parameter a which is chosen such that R! jr a ~ x remains 
finite in the limit of approach to the singularity along the trajectories coming out. In the 
a < 3 cases, we again get a family of RNGs coming out of the singularity along the apparent 
horizon. This happens because in general along the outgoing RNGs (here for simplicity we 
discuss only null geodesies) dR/dr 3 = (1 — \J f + F/R/ yj\ + f)R'/ (3r 2 ), so even with / ^ 
we again have situations with R = F + M s r 3+S , such that the first term in bracket goes to 
zero and the second term blows up but their product remains finite, giving the value of the 
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0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.00 0.01 0.02 0.03 0.04 0.05 

R R 

(a) (b) 



Figure 3: Different outgoing timelike geodesies for / = F = r 3 — 2.5r 4 . Density smoothly 
goes to zero at the boundary. All these geodesies show R = F kind of behaviour near the 
singularity. Different timelike geodesies meeting the same point correspond to the particle 
having different velocities, a) Different timelike singular geodesies meeting the apparent 
horizon at the boundary, b) Different timelike singular geodesies meeting R=0.02 at the 
boundary. 

tangent to be F . 

Finally, we make here some remarks in order to clarify the situation regarding the strength 
of the naked singularity for the cases we have considered here. As we indicated earlier, the 
singularity in this model is always gravitationally strong||. We discuss this in some detail 
below. 



The basic idea, in order to examine the strength of a singularity (naked or otherwise), is 
to examine the rate of curvature growth along the non-spacelike geodesies terminating in the 
singularity, in the limit of approach to the singularity. There are many criteria available for 
this purpose, as we point out below. According to the definitions of strength of a singularity 
as given by Tipler g |, |10| there are two criteria given to check the strength. According to 
one definition, a divergence condition to be satisfied is, RijV l V^ should go as 1/k 2 (k is the 
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affine parameter along the geodesic, vanishing in the limit of approach to the singularity) 
as one goes to the singularity, along every geodesic coming out, and also for the ingoing 
ones. According to the second definition of Tipler, the singularity is strong if this divergence 
condition is satisfied along at least one trajectory terminating in the singularity. We like to 
use this later definition of strength, because there is practically no way available to check the 
earlier definition operationally, because it is nearly impossible to integrate all non-spacelike 
geodesies in most of the collapse scenarios. Also, considering the complexity of the field 
equations, for most of the collapse models (consider e.g. the well-known case of Vaidya- 
Papapetrou radiation collapse), there is going to be some kind of a directional dependence 
always in the behaviour of the curvature growth in the limit of approach to the singularity. 
A uniform kind of curvature dependence in all directions does not seem plausible. Thus, 
it looks extremely reasonable to make the statement that the geodesic ends in a strong 
curvature singularity provided the curvatures grow as per the above requirement along it, 
rather than defining the strength as some kind of an absolute property of the singularity. 

In fact, there are other criteria as well proposed and available to check the strength. For 
example, according to the criterion given by Krolak (Tl| (see e.g. flOp, a divergence R i jV l V : > 
going as 1/k is sufficient to call the singularity strong. According to this criterion as well, 
clearly the naked singularities considered here are all strong curvature singularities along 
every null geodesic terminating in the singularity. There have also been some recent general 
considerations on strength by Nolan fll2| , according to which practically all shell- focusing 
singularities occurring in spherically symmetric spacetimes are gravitationally strong. Taking 
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all these considerations into account, it would appear quite reasonable to treat the naked 
singularities considered here to be strong curvature ones, which are physically important and 
not removable from the spacetime in a classical manner. 

4 Conclusion 

We have explored here the behaviour of null geodesies in the vicinity of the naked singularity 
developing in the gravitational collapse of a dust cloud. Together with the earlier information 
available in this connection (especially in the case n — 3), this provides a rather complete 
picture of how null trajectories traverse away from the naked singularity. 

We also see that there are ingoing and outgoing timelike geodesies terminating at the 
central naked singularity, and they have a behaviour very similar to that of the null geodesies 
when near the central singularity. This is an intriguing phenomenon and because of this in 
many cases it may become difficult to decide whether a given vector is null, timelike or space- 
like at the central singularity. This can happen because normally while doing calculations 
in such cases we finally use some approximations and take limits. 

We further showed that there is a family of radial null (as well as radial timelike geodesies) 
coming out of the singularity with the same ultimate tangent as that of the apparent horizon, 
when the parameter a < 3, i.e. the geodesies have R = F kind of behaviour as we approach 
the singularity. Very similar results can be obtained for spacelike geodesies also. For timelike 
geodesies in a = 3 case again we get results quite similar to that of null geodesies, as K l 
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has a power-law divergence near the singularity. Even in collapse scenarios more general 
than dust we can expect a similar phenomenon, i.e. in the cases corresponding to a < 3, 
we expect a family of geodesies coming out with R = F kind of behaviour as the general 
structure of dR/dr equation looks very similar. These results provide us with an insight into 
the causal structure in the vicinity of the naked singularity. 

We thank the referee for helpful comments on the paper. It is a pleasure to thank I. H. 
Dwivedi for discussions, and SSD would like to thank T. Harada, H. Iguchi, and K. Nakao 
for many useful comments. 
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